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The quantum master equation obtained by generalizing the geometric formulation of nonequi- 
librium thermodynamics to dissipative quantum systems is seriously nonlinear. We argue that 
nonlinearity occurs naturally in the step from reversible to irreversible equations and we analyze 
the nature of the nonlinear contribution. The thermodynamic nonlinearity leads to proper equilib- 
rium solutions and improves the dynamic behavior of dissipative quantum systems. We also discuss 
the Markovian character of the thermodynamic quantum master equation. The general ideas are 
illustrated for two-level systems. 
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Introduction. — For quantum information processing 
and a number of further present-day applications, dis- 
sipative quantum systems are a topic of key relevance. 
In a recent paper [1], thermodynamic arguments have 
been used to develop a master equation for quantum sys- 
tems that are weakly coupled to a classical environment. 
This quantum master equation obtained from thermody- 
namics is significantly different from the master equations 
usually formulated with the sole guidance from principles 
of quantum mechanics (2, 3]. For example, in the tradi- 
tional approach, the linear Lindblad equation is obtained 
under a Markovian assumption j4|. The thermodynamic 
quantum master equation is neither linear nor Marko- 
vian. It is hence important to elaborate and illustrate 
some of its general properties, and that is the purpose of 
the present letter. 

Whereas linearity seems natural in a quantum me- 
chanical setting, it is quite unusual in thermodynam- 
ics. This is a consequence of the appearance of en- 
tropy, which typically involves logarithmic terms. Go- 
ing beyond reversibility requires to go beyond linear- 
ity. As thermodynamics is the language for formulat- 
ing healthy equations with well-behaved solutions, non- 
linearity should not at all be considered as a drawback. 
Upon elaborating these general remarks in this letter, 
it should become increasingly clear that the thermody- 
namic quantum master equation is a most promising tool 
for all problems in quantum dissipation that have so far 
been treated by quantum master equations without a sys- 
tematic and complete thermodynamic background. The 
thermodynamic quantum master equation allows us to 
couple quantum subsystems to all kinds of environments, 
including time-dependent ones, and it comes with a cor- 
responding equation describing the influence of the quan- 
tum subsystem on the evolution of its environment. 

Thermodynamic quantum master equation. — Based on 
thermodynamic considerations [1], the following master 
equation for the evolution of the density matrix or statis- 
tical operator p on a Hilbert space % has been proposed 
to characterize a quantum subsystem in contact with an 
arbitrary classical nonequilibrium system acting as an 



environment: 
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(1) 



In this equation, h and fee are Planck's constant (divided 
by 2ir) and Boltzmann's constant, respectively. The first 
term describes the reversible contribution to the evolu- 
tion generated by the Hamiltonian H via the commuta- 
tor. All other terms are of irreversible nature and re- 
sult from a coupling to the environment of the quantum 
subsystem. They are expressed through commutators in- 
volving the self-adjoint coupling operators Qj so that the 
normalization condition, tr p = 1, is automatically pre- 
served in time. The subscript p on a quantum observable 
indicates the modified self-adjoint operator 
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which is basically the product of A and p, but with a 
compromise between placing p to the left or the right of 
A. 

Whereas the type of each coupling is given by the 
observable Qj, the strength of a coupling is expressed 
in a dissipative bracket [ , ] defined as a binary opera- 
tion on the space of observables for the classical environ- 
ment (throughout this letter, boldface bracket symbols 
are used to distinguish classical brackets from quantum 
commutators and anticommutators). If the equilibrium 
or nonequilibrium states of the environment are charac- 
terized by state variables x, classical observables are func- 
tions or functionals of x, and their evaluation at a partic- 
ular point of the state space is indicated by the subscript 
x. The classical observables H e and S e in Eq. ([T]) are the 
energy and the entropy of the environment, respectively. 
Dissipative brackets are bilinear, symmetric, [4,B C ] 3 = 
[B e , A e ] J , positive, [A c , Ac] 3 > 0, and satisfy the Leibniz 
or product rule, [A C B C , C e ] j = A C [B C , C C ] J + B e [A e , C e ] j , 
for arbitrary environmental variables A e , B e , and C c . 

Basic features. — In view of the definition ([2]) of A p , the 
second term in Eq. ([1} will, in general, be nonlinear in 
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p. This definition can be rewritten in a form that brings 
out the relationship to another possible compromise in 
placing p and extracts the nonlinearity, 



A p = \ ( A P + P A + A 'p) > A 'p 



(3) 

Setting A' p = corresponds to the linearization of the 
master equation ([T]) which has been proposed, but not 
recommended, in We hence may think of A' p as the 
origin of nonlinearity in the full thermodynamically con- 
sistent quantum master equation. The importance of A p 
stems from its occurrence in the canonical correlation 
(; } p (see Eq. (4.1.12) of fl), 



(A;B) ( 



[ tr(p x Ap 1 - x B)d\ = tr(A p B), 
Jo 



(4) 



which provides the key structural element in the ge- 
ometric formulation of irreversible dynamics for quan- 
tum systems [l|. The canonical correlation is symmetric, 
(A;B) = (B;A) , and positive, (A; A) > 0. Moreover, 
averages can be obtained as special cases of canonical 
correlations, (A) = tr(Ap) = tr(A p ) = (A; 1) . 

Whenever the state x of the environment varies in time, 
the strength of the coupling in the thermodynamic quan- 
tum master equation ([1]) becomes time-dependent. We 
then have a particular kind of non-Markovian behavior, 
as has been elaborated in the literature 0, 0] . It is ac- 
tually natural to expect that the environment changes 
in time because, in the thermodynamic approach [l|, the 
quantum master equation ([l} comes together with a cor- 
responding equation for the evolution of environmental 
observables, 
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In this equation, { , } and [ , ] are the Poisson and dissi- 
pative brackets of the classical system, respectively 043 ■ 
In addition to the properties listed above for the dissipa- 
tive brackets [ , ] J , energy conservation in the environ- 
ment (except for exchange of energy with the quantum 
subsystem) is guaranteed by the degeneracy requirement 
[j4 c , H c ] = for all observables A . The Poisson bracket 
is bilinear, antisymmetric, and satisfies the Leibniz rule 
as well as the Jacobi identity, where the latter expresses 
the time-structure invariance of the Poisson bracket. All 
these properties of classical Poisson brackets are also sat- 
isfied by their quantum counterparts, the commutators. 

If one looks only at the master equation (JTJ, the 
occurrence of time-dependent coefficients suggests non- 
Markovian behavior. If one looks at the coupled evolu- 
tion for the quantum subsystem and the classical environ- 
ment by the thermodynamically coupled set of Eqs. ([1]) 



and (J5J, however, the Markovian character of the de- 
scription is restored, provided that the total system is 
closed. By comparing Eqs. (JTJ and ((3J) one notes obvious 
exchange terms between the two subsystems. 

Equilibrium solution. — The nonlinearity of the thermo- 
dynamic quantum master equation is essential for the 
existence of a proper equilibrium solution. If, for an en- 
vironment in equilibrium with the temperature T c , we 
have the conditions 
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for all j, then we obtain the proper equilibrium solution 
to Eq. ©, 



p eq (X exp 



H 



(7) 



To obtain this result, we have made use of the identity 
[A, p] = [A pi In p] (see jlj) in the last term of the master 
equation (jj for A = Qj. 

The guaranteed existence of a proper equilibrium solu- 
tion is a major advantage of the thermodynamic quantum 
master equation. It is deeply linked to the nonlinearity 
of the master equation. 

Two-level system. — For a fc-state (or, fc-level) system, 
the underlying Hilbert space is a fc-dimensional complex 
vector space which, without loss of generality, we can take 
as C fc . The space of observables is the fc 2 -dimensional 
real vector space of self-adjoint k x fc-matrices with com- 
plex entries. For the two-level system, which has success- 
fully been used to describe both nuclear magnetic reso- 
nance and spontaneous emission in quantum optics (lflj . 
we choose the 2 x 2-unit matrix / and the three Pauli 
matrices 
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as the base vectors of the space of observables. More 
precisely, we express every self-adjoint complex 2x2- 
matrix A in terms of a real scalar a and a real three- 
vector a — (01,02,03), 

A = 0(a, a) = - (al + a 1 a 1 + a 2 a 2 + a 3 cr 3 ) . (9) 

Note that a is the trace of A. Commutators and anti- 
commutators can then conveniently be expressed as 

[A, B] = i0(O, a x b), (10) 

{A, B} — 0(af3 + a ■ b, /3a + ab). (11) 

From Eq. (jlOp . we obtain an identity for the frequently 
occurring double commutators, 



[A, [A. B}} = O(0, [« 2 1 - aa\ ■ b), 



(12) 



where 1 is the 3 x 3- unit matrix. From Eq. (|11[) , we 
obtain 



2tr(AB) =a/3 + a-b. 



(13) 
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Arbitrary functions / of an observable A can be calcu- 
lated with the formula 



f(A) = (/+ + /_,[/+-/_] a A 



with 



/+ = / 



a + a 



/-=/ 



(14) 
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where a = \a\. Equation ([Til can be verified by induction 
for powers of A and then generalized by Taylor expansion. 
From Eq. (fT4"f we further conclude that (a + a)/2 and 
(a — a)/2 must be the eigenvalues of A. 

As the density matrix has trace unity, it can be written 

as 



p=0(l,m) 



(16) 



For the eigenvalues to be nonnegative, we need m = 
\m\ < 1. This set of admissible choices of m is known as 
the Bloch sphere. For m = 1, one of the two eigenvalues 
of p is zero and we have a pure state. From Eq. fllUl) . 
we obtain {A) p = (a + a ■ m)/2, which implies that the 
jth component of m is given by the average (<7j) . By 



using Eqs. (ITUl) and (fT4")l in Eq. ©, we find the following 
explicit form for the nonlinear part of A p , 



A' = -O(0,p(m)[m 2 1 - mm] ■ a 



with 



p(m) 



1 



m- 
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(17) 



(18) 



The function p,(m) is displayed in Figure [JJ The singu- 
larities of the two terms in Eq. (|18l) at m = cancel so 
that p(m) « 1/3 for small m. According to Eq. (|12[) . the 



factor [m 1 — mm] may be regarded as a double commu- 
tator formed with p. The nonlinear contribution to the 
quantum master equation associated with p(m) leads to 
an improved relaxation behavior, as we shall see below. 

We now choose the Hamiltonian H = O(0,hujq 3 ), 
where uj is the angular frequency associated with the 
energy difference between the two levels of the system 
and q 3 = (0,0,1), as well as the two coupling operators 
Qj = O(0, qj) with q l = (1,0,0) and q 2 = (0,1,0). The 
environment be a heat bath, the state of which be char- 
acterized by its energy H e . The thermodynamic prop- 
erties of the bath are given by the thermodynamic re- 
lationship S e (H e ), and we introduce the bath tempera- 
ture by \/T c = dS c /dH c . This temperature characterizes 
the black-body radiation to which our system is exposed. 
Both dissipative brackets are assumed to be given by 



[A c ,B e ] J = _ 7o ___ 



dH, 
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(19) 




FIG. 1: The function p(m) characterizing the nonlinear con- 
tribution to the thermodynamic quantum master equation for 
a two- level system [see Eq. (|18[) ]. 



The quantum master equation {JJ can now be recog- 
nized to be equivalent to an evolution equation for m, 
known as the Bloch equation (Tl| . 



dm 
~dt 
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- ^foq3+lo^{m 2 l + mm)-q 3 , (20) 

with R = (1 + q 3 q 3 )/2. Our choice of the coupling op- 
erators Qj is motivated by the two Lindblad operators 
that have been derived for quantum optical applications 
of two-level systems for the case of spontaneous emission 
(see, for example, Eq. (3.219) of 0). As an alternative, 
one could include Q 3 — O(0, q 3 ) as a third operator with 
the same dissipative bracket (|19[) to achieve an isotropic 
frictional coupling. The only effect would be to change 
the anisotropic matrix R in Eq. ([2"U| into the unit ma- 
trix 1 . This would correspond to the case of strong colli- 
sions that cause energy decay whenever they cause dipole 
phase interruption [12j . For nuclear spin relaxation, for 
which the Bloch equation had been proposed originally, 
the situation with R — 1 can also be realized, namely in 
isotropic molecular environments, both in gases and in 
low- viscosity liquids [13]. It is well known, however, that 
longitudinal relaxation rates that are by orders of mag- 
nitude smaller than the transverse ones are much more 
typical for nuclear spin relaxation This situation 

can be achieved by enhancing the coupling strength as- 
sociated with q 3 to become the dominating one. 
The equilibrium solution of Eq. (|20p is given by 



-q 3 tanh 



froj 
2k B T c 



(21) 



where 70 is the spontaneous emission rate. 



Contrary to the steady solution of the linearized quan- 
tum master equation, m = —q^fkjj/(2k B ,T^), the solution 
([2T]) always lies in the Bloch sphere, even for very low 
temperatures. The time-dependent solution of the ther- 
modynamic quantum master equation can actually never 
leave the Bloch sphere, which is a nice consequence of 
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thermodynamic consistency (note that, for very low tem- 
peratures, the solution of the linearized master equation 
leaves the Bloch sphere; although such low temperatures 
are unrealistic, this illustrates the improved relaxation 
behavior due to the nonlinearity) . For a small deviation 
m' from the steady state solution ([2T]) . we obtain a lin- 
earized version of Eq. (|20[) , 



dm 1 
~dT 



ujq 3 x m 
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70 — (1 + 3q 3 g3j • m - 70m — q 3 <?3 ■ m 



2k B T e , 
70 — R ■ m 



dm 



(22) 



where p{m) and its derivative are to be evaluated at m eq . 
The nonlinear terms enhance the relaxation, most dra- 
matically near the boundary of the Bloch sphere. 

The nonlinear part of A p in Eq. (| 1 T[) can be rewritten 
in the alternative form 



A' = /j,(m) Apti(AAp) - Ati(Apf 



(23) 



where Ap = p — 1/2 is the deviation of the density ma- 
trix from a uniform distribution over all states (which 
corresponds to the entropy-dominated high-temperature 
limit). This reformulation suggests that the nonlinear- 
ity favors a uniform distribution and that the strength of 
this effect is proportional to p(m). This is an effect of 
quantum fluctuations; it is independent of temperature 
and disappears for classical systems. 

We had argued before that, in view of the presence 
of entropy, nonlinearity occurs very naturally in thermo- 
dynamics. How can the nonlinearity then disappear in 
the classical limit? This is indeed a fortuitous cancela- 
tion occurring for classical systems on the level of dis- 
tribution functions because of the identity fd(5S/Sf) = 



— fce/rfln/ = —k^df (it is shown in Section III.B of Q 
that, in a thermodynamic approach, this property leads 
to the linearity of the Fokker-Planck equation for classi- 
cal systems) . The noncommutativity of quantum observ- 
ables prevents such a cancelation. 

Summary and outlook. — In this letter, we have estab- 
lished some advantages of the thermodynamically consis- 
tent quantum master equation. The resulting equation 
is significantly different from the ones commonly con- 
structed without guidance from thermodynamics. In par- 
ticular, the thermodynamic master equation is nonlinear 
and this nonlinearity is very helpful for obtaining physi- 
cally meaningful solutions (for example, at equilibrium). 

The thermodynamic quantum master equation de- 
scribes the influence of any classical environment on a 
quantum subsystem. Moreover, it is supplemented by an 
equation describing the reverse influence of the quantum 
subsystem on the environment. If the total system is 
closed, we obtain a Markovian description of the coupled 
subsystems even if the coefficients in the quantum master 
equation change with a changing environment. 

As quantum master equations are nowadays employed 
in many applications involving dissipative quantum sys- 
tems, the nonlinear thermodynamic quantum master 
equation offers a new perspective on many problems. 
Problems that involve more complicated environments 
than simple heat baths can be approached in a thermo- 
dynamically consistent way. 

Our thermodynamic analysis is currently restricted to 
quantum master equations. Also for other approaches to 
dissipative quantum systems, such as operator Langevin 
equations, stochastic dynamics in Hilbert space, or path 
integrals (3.I3. [T3 - Il^ . thermodynamic consistency should 
be established — this is work in progress. 
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